Observation of r) c decay into S+S and S H+ final states 

M. Ablikim 1 , M. N. Achasov 5 , O. Albayrak 3 , D. J. Ambrose 39 , F. F. An 1 , Q. An 40 , J. Z. Bai 1 , Y. Ban 27 , J. Becker 2 , 
J. V. Bennett 17 , M. Bertani 18A , J. M. Bian 38 , E. Boger 20 ' a , 0. Bondarenko 21 , I. Boyko 20 , R. A. Briere 3 , V. Bytev 20 , X. Cai 1 , 
O. Cakir 35A , A. Calcaterra 18A , G. F. Cao 1 , S. A. Cetin 35S , J. F. Chang 1 , G. Chelkov 20 '", G. Chen 1 , H. S. Chen 1 , 
J. C. Chen 1 , M. L. Chen 1 , S. J. Chen 25 , X. Chen 27 , Y. B. Chen 1 , H. P. Cheng 14 , Y. P. Chu 1 , F. Coccetti 18A , 
D. Cronin-Hennessy 38 , H. L. Dai 1 , J. P. Dai 1 , D. Dedovich 20 , Z. Y. Deng 1 , A. Denig 19 , I. Denysenko 20 ' 6 , M. Destefanis 43A ' 43C , 
W. M. Ding 29 , Y. Ding 23 , L. Y. Dong 1 , M. Y. Dong 1 , S. X. Du 46 , J. Fang 1 , S. S. Fang 1 , L. Fava 43S ' 43C , F. Feldbauer 2 , 
C. Q. Feng 40 , R. B. Ferroli 18A , C. D. Fu 1 , J. L. Fu 25 , Y. Cao 34 , C. Ceng 40 , K. Goetzen 7 , W. X. Gong 1 , W. Gradl 19 , 
M. Greco 43A - 43C , M. H. Gu\ Y. T. Gu 9 , Y. H. Guan 6 , A. Q. Guo 26 , L. B. Guo 24 , Y. P. Guo 26 , Y. L. Han 1 , F. A. Harris 37 , 
K. L. He 1 , M. He 1 , Z. Y. He 26 , T. Held 2 , Y. K. Heng 1 , Z. L. Hou 1 , H. M. Hu\ J. F. Hu 36 , T. Hu\ G. M. Huang 15 , 

G. S. Huang 40 , J. S. Huang 12 , X. T. Huang 29 , Y. P. Huang 1 , T. Hussain 42 , C. S. Ji 40 , Q. Ji\ Q. P. Ji 26,c , X. B. Ji\ 
X. L. Ji 1 , L. L. Jiang 1 , X. S. Jiang 1 , J. B. Jiao 29 , Z. Jiao 14 , D. P. Jin 1 , S. Jin 1 , F. F. Jing 34 , N. Kalantar-Nayestanaki 21 , 
M. Kavatsyuk 21 , M. Kornicer 37 , W. Kuehn 36 , W. Lai 1 , J. S. Lange 36 , C. H. Li 1 , Cheng Li 40 , Cui Li 40 , D. M. Li 46 , F. Li 1 , 
G. Li 1 , H. B. Li 1 , J. C. Li 1 , K. Li 10 , Lei Li 1 , Q. J. Li 1 , S. L. Li 1 , W. D. Li 1 , W. G. Li 1 , X. L. Li 29 , X. N. Li 1 , X. Q. Li 26 , 

X. R. Li 28 , Z. B. Li 33 , H. Liang 40 , Y. F. Liang 31 , Y. T. Liang 36 , G. R. Liao 34 , X. T. Liao 1 , B. J. Liu 1 , C. L. Liu 3 , C. X. Liu 1 , 

C. Y. Liu 1 , F. H. Liu 30 , Fang Liu 1 , Feng Liu 15 , H. Liu 1 , H. H. Liu 13 , H. M. Liu 1 , H. W. Liu 1 , J. P. Liu 44 , K. Y. Liu 23 , 
Kai Liu 6 , P. L. Liu 29 , Q. Liu 6 , S. B. Liu 40 , X. Liu 22 , Y. B. Liu 26 , Z. A. Liu 1 , Zhiqiang Liu 1 , Zhiqing Liu 1 , H. Loehner 21 , 

G. R. Lu 12 , H. J. Lu 14 , J. G. Lu 1 , Q. W. Lu 30 , X. R. Lu 6 , Y. P. Lu\ C. L. Luo 24 , M. X. Luo 45 , T. Luo 37 , X. L. Luo 1 , M. Lv\ 

C. L. Ma 6 , F. C. Ma 23 , H. L. Ma 1 , Q. M. Ma 1 , S. Ma 1 , T. Ma 1 , X. Y. Ma 1 , Y. Ma 11 , F. E. Maas 11 , M. Maggiora 43A ' 43C , 
Q. A. Malik 42 , Y. J. Mao 27 , Z. P. Mao 1 , J. G. Messchendorp 21 , J. Min 1 , T. J. Min 1 , R. E. Mitchell 17 , X. H. Mo 1 , C. Morales 

Morales 11 , C. Motzko 2 , N. Yu. Muchnoi 5 , H. Muramatsu 39 , Y. Nefedov 20 , C. Nicholson 6 , I. B. Nikolaev 5 , Z. Ning 1 , 
S. L. Olsen 28 , Q. Ouyang 1 , S. Pacetti 18S , J. W. Park 28 , M. Pelizaeus 2 , H. P. Peng 40 , K. Peters 7 , J. L. Ping 24 , R. G. Ping 1 , 
R. Poling 38 , E. Prencipe 19 , M. Qi 25 , S. Qian 1 , C. F. Qiao 6 , X. S. Qin 1 , Y. Qin 27 , Z. H. Qin 1 , J. F. Qiu 1 , K. H. Rashid 42 , 
G. Rong 1 , X. D. Ruan 9 , A. Sarantsev 20 ' d , B. D. Schaefer 17 , J. Schulze 2 , M. Shao 40 , C. P. Shen 37 ' e , X. Y. Shen 1 , 

H. Y. Sheng 1 , M. R. Shepherd 17 , X. Y. Song 1 , S. Spataro 43A ' 43C , B. Spruck 36 , D. H. Sun 1 , G. X. Sun 1 , J. F. Sun 12 , 
S. S. Sun 1 , Y. J. Sun 40 , Y. Z. Sun 1 , Z. J. Sun 1 , Z. T. Sun 40 , C. J. Tang 31 , X. Tang 1 , I. Tapan 35C , E. H. Thorndike 39 , 

D. Toth 38 , M. Ullrich 36 , G. S. Varner 37 , B. Wang 9 , B. Q. Wang 27 , D. Wang 27 , D. Y. Wang 27 , K. Wang 1 , L. L. Wang 1 , 
L. S. Wang 1 , M. Wang 29 , P. Wang 1 , P. L. Wang 1 , Q. Wang 1 , Q. J. Wang 1 , S. G. Wang 27 , X. L. Wang 40 , Y. D. Wang 40 , 

Y. F. Wang 1 , Y. Q. Wang 29 , Z. Wang 1 , Z. G. Wang 1 , Z. Y. Wang 1 , D. H. Wei 8 , J. B. Wei 27 , P. Weidenkaff 19 , Q. G. Wen 40 , 

S. P. Wen 1 , M. Werner 36 , U. Wiedner 2 , L. H. Wu\ N. Wu 1 , S. X. Wu 40 , W. Wu 26 , Z. Wu\ L. G. Xia 34 , Z. J. Xiao 24 , 
Y. G. Xie 1 , Q. L. Xiu\ G. F. Xu\ G. M. Xu 27 , H. Xu\ Q. J. Xu 10 , X. P. Xu 32 , Z. R. Xu 40 , F. Xue 15 , Z. Xue 1 , L. Yan 40 , 



W. B. Yan 40 , Y. H. Yan 16 , H. X. Yang 1 , Y. Yang 15 , Y. X. Yang 8 , H. Ye 1 , M. Ye 1 , M. H. Ye 4 , B. X. Yu 1 , C. X. Yu 26 , 

H. W. Yu 27 , J. S. Yu 22 , S. P. Yu 29 , C. Z. Yuan 1 , Y. Yuan 1 , A. A. Zafar 42 , A. Zallo 18A , Y. Zeng 16 , B. X. Zhang 1 , 
B. Y. Zhang 1 , C. Zhang 25 , C. C. Zhang 1 , D. H. Zhang 1 , H. H. Zhang 33 , H. Y. Zhang 1 , J. Q. Zhang 1 , J. W. Zhang 1 , 
J. Y. Zhang 1 , J. Z. Zhang 1 , S. H. Zhang 1 , X. J. Zhang 1 , X. Y. Zhang 29 , Y. Zhang 1 , Y. H. Zhang 1 , Y. S. Zhang 9 , 
Z. P. Zhang 40 , Z. Y. Zhang 44 , G. Zhao 1 , H. S. Zhao 1 , J. W. Zhao 1 , K. X. Zhao 24 , Lei Zhao 40 , Ling Zhao 1 , M. G. Zhao 26 , 
Q. Zhao 1 , Q. Z. Zhao 9 ^, S. J. Zhao 46 , T. C. Zhao 1 , X. H. Zhao 25 , Y. B. Zhao 1 , Z. G. Zhao 40 , A. Zhemchugov 20 '", B. Zheng 41 , 
J. P. Zheng 1 , Y. H. Zheng 6 , B. Zhong 24 , J. Zhong 2 , Z. Zhong 9J ', L. Zhou 1 , X. K. Zhou 6 , X. R. Zhou 40 , C. Zhu 1 , K. Zhu\ 
K. J. Zhu 1 , S. H. Zhu 1 , X. L. Zhu 34 , Y. C. Zhu 40 , Y. M. Zhu 26 , Y. S. Zhu 1 , Z. A. Zhu 1 , J. Zhuang 1 , B. S. Zou 1 , J. H. Zou 1 

(BESIII Collaboration) 

1 Institute of High Energy Physics, Beijing 10004-9, P. R. China 
2 Bochum Ruhr- University, 44780 Bochum, Germany 
3 Carnegie Mellon University, Pittsburgh, PA 15213, USA 
4 China Center of Advanced Science and Technology, Beijing 100190, P. R. China 
5 G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 

6 Graduate University of Chinese Academy of Sciences, Beijing 100049, P. R. China 

7 GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 

8 Guangxi Normal University, Guilin 541004, P. R. China 
9 GuangXi University, Naming 530004, P.R.China 
10 Hangzhou Normal University, Hangzhou 310036, P. R. China 
11 Helmholtz Institute Mainz, J.J. Becherweg 45, D 55099 Mainz,Germany 
12 Henan Normal University, Xinxiang 453007, P. R. China 
13 Henan University of Science and Technology, Luoyang 471003, P. R. China 
14 Huangshan College, Huangshan 245000, P. R. China 
15 Huazhong Normal University, Wuhan 430079, P. R. China 
16 Hunan University, Changsha 410082, P. R. China 
17 Indiana University, Bloomington, Indiana 47405, USA 
18 (A)INFN Laboratori Nazionali di Frascati, Frascati, Italy; (B)INFN and University of Perugia, 1-06100, Perugia, Italy 
19 Johannes Gutenberg University of Mainz, Johann-Joachim-Becher-Weg 45, 55099 Mainz, Germany 
20 Joint Institute for Nuclear Research, 141980 Dubna, Russia 
21 KV1 '/University of Groningen, 9747 AA Groningen, The Netherlands 
22 Lanzhou University, Lanzhou 730000, P. R. China 



23 Liaoning University, Shenyang 110036, P. R. China 
24 Nanjing Normal University, Nanjing 210046, P- R- China 

25 Nanjing University, Nanjing 210093, P. R. China 

26 Nankai University, Tianjin 300071, P. R. China 

27 Peking University, Beijing 100871, P. R. China 

28 Seoul National University, Seoul, 151-747 Korea 

29 Shandong University, Jinan 250100, P. R. China 

30 Shanxi University, Taiyuan 030006, P. R. China 
31 Sichuan University, Chengdu 610064, P. R- China 

32 Soochow University, Suzhou 215006, China 
33 Sun Yat-Sen University, Guangzhou 510275, P. R. China 
34 Tsinghua University, Beijing 100084, P- R- China 
(A) Ankara University, Ankara, Turkey; (B)Dogus University, Istanbul, Turkey; (C)Uludag University, Bursa, Turkey 

36 Universitaet Giessen, 35392 Giessen, Germany 
37 University of Hawaii, Honolulu, Hawaii 96822, USA 
38 University of Minnesota, Minneapolis, MN 55455, USA 
39 University of Rochester, Rochester, New York 14627, USA 
40 University of Science and Technology of China, Hefei 230026, P. R. China 
41 University of South China, Hengyang 421001, P. R. China 
42 University of the Punjab, Lahore-54590, Pakistan 
(A)University of Turin, Turin, Italy; (B)University of Eastern Piedmont, Alessandria, Italy; (C)INFN, Turin, Italy 

44 Wuhan University, Wuhan 430072, P. R. China 
45 Zhejiang University, Hangzhou 310027, P. R. China 
46 Zhengzhou University, Zhengzhou 450001, P. R. China 

a also at the Moscow Institute of Physics and Technology, Moscow, Russia 
b on leave from the Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 
c Nankai University, Tianjin, 300071, China 
d also at the PNPI, Gatchina, Russia 
e now at Nagoya University, Nagoya, Japan 
1 Guangxi University, Nanning, 530004, China 



Using a data sample of 2.25 x 10 8 J ftp events collected with the BESIII detector, we present the first 
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observation of the decays of r/ c mesons to E E and 3 3 . The branching fractions are measured 
to be (2.11±0.28 st at. ±0.18 syst . ±0.50pdg) x 10" 3 and (0.89±0.16 sta t. ±0.08 syst . ±0.21 PDG ) x HT 3 
for 77 c — > E + E~ and H~H + , respectively. These branching fractions provide important information 
on the helicity selection rule in charmonium-decay processes. 

PACS numbers: 13.25.Gv, 13.20.Gd, 14.40.Pq 



I. INTRODUCTION 

Experimental studies on exclusive charmonium decays play an important role in testing perturbative Quantum 
Chromodynamics (pQCD). In the Standard Model (SM), the rj c meson is the lowest lying charmonium state in a 
h spin-parity configuration. Although the n c cannot be produced directly from e + e~ annihilations, it is produced 
copiously in radiative decays of J/ip and tp' [l|. The large Jji\> and i\>' data samples taken with the BESIII detector 
at the BEPCII provide an opportunity for a detailed study of r\ c decays. 

The complexity of QCD remains unsolved in the charmonium-mass region, and there are still many contradictions 
between pQCD calculations and experimental measurements. In particular, the pQCD helicity selection rule 2j-|4| is 
violated in many exclusive charmonium-decay processes, for example, the decay processes with meson pairs in the 
final state, like J/ip — )• VP, rj c —> VV, and Xd —> VV, where V and P denote vector and pseudoscalar mesons. 
Other examples include decay processes with baryon anti-baryon pairs in the final state, such as r\ c — > BgBg, and 
Xco ~^ BgBg, where BgBg denote the octet baryon anti-baryon pairs. Many attempts have been made to understand 
these contradictions, such as by the quark-diquark model for the proton Q,Q|, constituent quark-mass corrections 7, 8|, 
mixing between the charmonium state and the glueball [9J, and the quark pair creation model [101 ] - However, the 
measured branching fractions are not consistent with the predictions of any of these models. 

In Refs. 
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12j . intermediate meson loop (IML) transitions are proposed, where the long-distance interaction 
can evade the Okubo-Zweig-Iizuka (OZI) rule and allow the violation of the pQCD helicity selection rule. Further 
calculations on the branching fractions of r\ c — \ B^B^, Xco Bg,Bg and h c — > BgB^ based on charmed-meson loops 



were carried out 13J, and the results agree with the measured branching fractions of rj c — > pp and rj c — > A A. Using a 
sample of 2.25 x 10 s J/tp events |l4[ collected with the BESIII detector in 2009, we measure the branching fractions 
of T] c — > E + S _ and rj c — > for the first time via the J/ip —> ^r\ c radiative decay process. 



■5 



II. DETECTOR AND MONTE CARLO SIMULATION 



BEPCII 15| is a double-ring e + e collider designed to provide a peak luminosity of 10 33 cm 2 s 1 at a center-of- 



mass energy of 3.77 GeV. The BESIII 15| detector has a geometrical acceptance of 93% of 47r and has four main 
components: (1) A small-cell, helium-based (40% He, 60% C3H8) main drift chamber (MDC) with 43 layers providing 
an average single- hit resolution of 135 /im, charged-particle momentum resolution in a 1 T magnetic field of 0.5% 
at 1 GeV/c. (2) An electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) crystals in a cylindrical structure 
(barrel) and two endcaps. For 1 GeV photons, the energy resolution is 2.5% (5%) in the barrel (endcaps), and the 
position resolution is 6 mm (9 mm) in the barrel (endcaps). (3) A time-of-flight system (TOF) consisting of 5-cm-thick 
plastic scintillators, with 176 detectors of 2.4 m length in two layers in the barrel and 96 fan-shaped detectors in the 
endcaps. The barrel (endcaps) time resolution of 80 ps (110 ps) provides 2er K/tt separation for momenta up to 
~ 1 GeV/c. (4) The muon system (MUC) consists of 1000 m 2 of Resistive Plate Chambers (RPCs) in 9 barrel and 8 
endcap layers and provides 2 cm position resolution. 

The optimization of the event selection and the estimate of backgrounds are performed using Monte Carlo (MC) 
simulated data. The GEANT4 [16j]-based simulation software BOOST ljj includes the geometry and the material 
description of the BESIII spectrometer, the detector response and digitization models, as well as the tracking of the 
detector running conditions and performances. The production of the J / tj: resonance is simulated by the MC event 



generator KKMC 
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19( , while the decays are generated by EvtGen 



201 ] for the known decay modes with branching 



fractions set to world average values [lj, and by LundCharm [2l| for the remaining unknown decays. 



III. EVENT SELECTION 



We select rj c mesons via the radiative decay J/t/j — > "/rj c with its subsequent decay into and The 

S + candidates are reconstructed from the decay S + — > pn° with the 7r° decaying into a pair of photons; the S - 
candidates are reconstructed from the decays S~ — » Att~ and A — > pir~ . The anti-particle candidates, S~ and S + , 
are reconstructed in a similar way but with the decay products changed to the corresponding anti-particles. 

Tracks of charged particles in the polar-angle range |cos#| < 0.93 are reconstructed from hits in the MDC. The 
TOF and dE/dx information are combined to form particle identification (PID) confidence levels for the 7T, K and p 
hypotheses. Each track is assigned to the particle type that corresponds to the hypothesis with the highest confidence 
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level. Photon candidates are reconstructed by clustering the energy deposited in the EMC crystals. The minimum 
energy requirement is 25 MeV for barrel showers (| cos#| < 0.80) and 50 MeV for endcap showers (0.86 < | cos#| < 
0.92). Requirements on the EMC cluster timing are applied to suppress electronic noise and energy deposits unrelated 
to the event. Candidate 7r° mesons are reconstructed from pairs of photons with an invariant mass in the range 
0.115 GeV/c 2 < M(77) < 0.155 GeV/c 2 . The tt° invariant-mass resolution is determined to be 4.2 MeV/c 2 by fitting 
the invariant-mass distribution of the 77 pairs from data after applying all the requirements except for the 7r°-mass 
window, as shown in Fig.[lja). In the fit, the ir° signal is taken with a Gaussian form, and the background is described 
by a second-order Chebychev polynomial function. 




FIG. 1: (a) A fit to the invariant-mass distribution of 77 pairs after applying all the requirements except for the 7r°-mass 
window. Dots with error bars are data, and the solid line is the total fit result. The signal is represented by the short-dashed 
line and the background by the long-dashed line, (b) A scatter plot for M(p7r°) versus M(p7r°). (c) Invariant-mass distributions 
of pir° and pn°; solid dots with error bars are M(p-7r°), and the open circles with error bars are M(p7r°). 
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FIG. 2: (a) Scatter plot for M(p7r + 7r + ) versus M(p7r n ). Invariant-mass distributions of (b) pn and pn + , and (c) pir n and 
p7r + 7r + . Solid dots with error bars are M(p7r~) and M(p7r~7r~), and open circles with error bars are M(p7r + ) and M(piv + -K + ). 



For J/ip — > j7] c — > 7^+E , exactly one proton, one anti-proton, at least five photons and at least two ir° candidates 
from the combination of these photons are required. A four-constraint (4C) kinematic fit, based on momentum and 
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energy conservation, is applied under the J/ip — > 7pp7r°7r° hypothesis, and \\c < 30 is required. For events with more 
than five photons or more than two ir° candidates, the combination with the minimum x\c * s retained in the analysis. 
The events are also fitted to the J/tp — > ppir°Tr and J/ip — > r yjppTr°ir° hypotheses. We require X4c(PP 7r ° 7r °) > ^00 
and X4C"(7Pp7r°7r°) < X4c(77PP 7r ° 7r °)- The p, p and the two ir° candidates are combined to form the E + and S - 
candidates by minimizing (M p7r o— M s +) 2 + (M^o — M^-) 2 . Furthermore, the combined p, tt° (p, tt°) pair must have 
an invariant mass within 15 MeV/c 2 of the S + mass, as shown in Fig. [ljb) and (c). 

For J/tp — >• jr/c 7S _ S + , exactly one proton, one anti-proton, two 7r + s, two 7r _ s and at least one photon are 
required. A 4C kinematic fit is applied under the J/tp —t r yppK + ir + ir~ ir~ hypothesis, and x\c < 90 is required. 
For events with more than one photon candidate, only the combination with the minimum x\c ^ s retained in the 
analysis. The events are also fitted to the J/tp — > pp7r + 7r + 7r~7r~ and J/tp — > r y"/ppTt + n + n~ n~ hypotheses. We require 
xlc(PP 7r+7r+7r ^ 7r ^) > 200 and X4c(7M" r+7r+7r ~ Tr ~) < X4c(TTPP 7r+7r+7r_7r_ )- 

To reconstruct the kinematical information of A and S _ , vertex fits are applied to the charged tracks (jyir~ and 
p7r~7r~ for A and S - , respectively), with the requirement that all the tracks originated from the same decay point. 
Next, secondary vertex fits are applied to these reconstructed particles, with the requirement that their flight time is 
consistent with the one predicted from their final-state particles. The p, n~ (p, n + ) combination with an invariant 
mass that is the closest to the A (A) mass is chosen to form the A (A). Furthermore, the mass difference must be 
within 10 MeV/c 2 , as shown in Fig.(2jb). The p, ir~ , ir~ (p, 7r + , ir + ) combination must have an invariant mass within 
9 MeV/c 2 of the S~ mass, as shown in Fig.[2ja) and (c). 

Figure |3] shows the invariant-mass distributions of and S _ S + pairs after applying all the event selection 

criteria. A clear signature of an i] c resonance is observed. 

IV. BACKGROUND STUDIES 

The background can be classified into two categories: background from rj c decays which produces a peak within the 
rj c signal region, and background from J/tp decays which gives a smooth distribution under the rj c resonance. 

For r\ c — > the potential peaking background channel is ij c — > ppir°ir , which has not previously been measured. 

By requiring the invariant mass of any pir° combination to be outside a mass window of 50 MeV/c 2 centered at the 
S + mass and the ppir°ir° invariant mass within 30 MeV/c 2 from the r\ c mass, the number of rj c — > ppir°ir events is 
obtained, and the branching fraction is determined to be (5.0 ± 0.6 s t a t.) x 10 -4 , where the uncertainty is statistical 
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only. Out of 5 x 10 5 J/tp — > ^r\ c — > ^pp7T°TT MC simulated events, 193 events survive after applying the event selection 
criteria. Using the measured branching fraction, the background contribution from this process is estimated to be 
0.7 events. For the background from J/tp decays, the main sources are J/tp — > S + E~ and J/tp — > 7r°5] + S - , which 
have a fake photon or a photon from ir° that escaped from detection, respectively; and J/tp — > 7S + S _ , which is 
an irreducible background to the signal process. Using 5 x 10 5 MC simulated events for each channel and applying 
the event selection criteria to these MC samples, the background contributions are estimated by normalizing the 
number of the surviving events to the total number of J/tp events. In the normalization, the branching fraction 
of J/tp — > is taken from Ref. [lj, and the branching fractions of J/tp — > 7£ + S~ and J/tp — > 7r°S + E _ are 

measured in this analysis. The branching fraction of J/tp — > 7r°£ + £~ is measured to be (5.0 ± 0.1 s t a t.) X 10~ 4 using 
similar event selection criteria but with an additional photon and a tt° reconstructed from the selected photons. The 
branching fraction of J/tp — >• 7£ + E~ is measured to be (7.4 ± 0.6 s t a t.) x 10 -5 with the same event selection criteria 
as was applied for the signal events, but without requiring that the system forms an rj c resonance and with 

a selection on the invariant mass of 2.4 GeV/c 2 < M(S + E~) < 2.8 GeV/c 2 . The total background is estimated to 
be 351 events in the entire mass region, as shown in Fig. |3Ia). The total background shape is found to be smooth 
without an enhancement under the rj c resonance. 




FIG. 3: Invariant mass distributions of data and MC background channels together with the fitted curves for (a) £+£", and 
(b) H~H + . Dots with error bars are data, and the histograms are the backgrounds from simulated J/tp decays. Solid lines 
are the total fit results, signals are shown in short-dashed lines, and backgrounds are shown as long-dashed lines and shaded 
histograms. 

For rj c — > the potential peaking background channels are rj c — ¥ ppTT + tt + tt~ 7r _ and rj c hkn + iT~ . Out of 
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2.5 x 10 5 simulated MC events for each channel, 2 and 21 events survived after applying the event selection criteria. 
The branching fractions of these two channels are determined to be (6.7 ± 1.0 s tat.) x 10 -4 and (6.3 ± 0.4 sta t.) x 10 -3 , 
respectively, where the uncertainties are statistical only. The invariant-mass requirements for rj c — ¥ ppir + Tr + ir~Tr~ 
are: |M p7r - — Ma| > 20 MeV/c 2 (no pn~ combination consistent with a A), IM^-^-— Ms- | > 25 MeV/c 2 (no p7r~7r~ 
combination consistent with a S - ), and |M J) p 7r + 7r + 7r - 7r - — M^ c | < 30 MeV/c 2 ; for r/ c — > AAit + it~, the only change 
is |M p7r - — Ma | < 20 MeV/c 2 . Using the measured branching fractions, the background contributions from the two 
peaking background channels are estimated to be 0.02 and 2 events to the signal after normalizing the number of the 
surviving events to the total number of the J/tp events, respectively. The main background channels from J/tp decays 
are J/tp — > S~S + and J/tp — > 7r°S - S + , which have one fake photon or one photon from the tt° that escaped from 
detection, and J/tp — > 7S~S + , which is an irreducible background to the signal. Another background contribution 
from J/tp — ► £°A7r + 7r~ — > jAAtt + tt^ + c.c. is apparently seen from the invariant-mass distribution of 7A pairs. To 
estimate the background contribution from the process J/tp — > 7r°S~S + including intermediate states, J/tp — > 7r°S - S + 
decays are reconstructed from data, and the signal yield is obtained in each M(S _ S + ) mass bin. The selection criteria 
are similar to that for signal events but with an additional photon and a ir° reconstructed from the selected photons. 
The relative efficiencies of the 75~S + and 7r°S - S + selection criteria are estimated in each M(5~5 + ) mass bin using 
J/tp — >■ 7r°S~S + MC events. Combining this relative efficiency with the number of J/tp — > 7r°S"S + signal events in 
each M(S~S + ) mass bin, the number of 7r S _ E! + events that pass the 7S~S + selection is estimated. We generated 
5 x 10 6 MC events for the channels J/ip S~S+ and J/tp -> S°A7r + Tr _ + c.c. and 2.5 x 10 5 MC events for the 
channel of J/tp — > 7S _ S + , and applied the event selection criteria to these MC samples. The contribution from each 
background process is estimated by normalizing the number of the surviving events to the total number of the J/tp 
events. In the normalization, the branching fraction of J/tp — > is taken from Ref. [l| and the branching fractions 

of J/tp — > 7S~S + and J/tp —> E°A7r + 7r~ are measured in this analysis. The branching fraction of J/tp — > £°A7r + 7r~ 
is measured to be (4.7 ± 0.1 s t a t.) x 10~ 4 by fitting the invariant-mass distribution of 7A pairs. The branching fraction 
of J/tp — > 7S~S + is measured to be (1.8 ± 0.5 s t a t.) x 10~ 5 by excluding the S~S + system to form an r\ c meson via 
the requirement M(S~S + ) < 2.8 GeV/c 2 . The total background from J/tp decays is estimated to be 116 events in the 
entire mass region, as shown in Fig. |3|h) , and is smoothly distributed and no enhancement under the r\ c resonance is 
observed. 
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V. SIGNAL EXTRACTIONS AND BRANCHING FRACTION CALCULATIONS 

Signal yields are obtained from unbinned maximum likelihood fits to the invariant-mass distributions of E + £~ and 
H _ !E! + candidates. The probability density function (PDF) used in the fit is given by 

F(m) = <r res ® (e(m) x Et x damping{E 1 ) x BW{m)) + BKG(m), 

where BW(m) and BKG(m) are the signal component described by the Breit-Wigner form and the background 
component, respectively; a res is the experimental resolution function and e(m) is the mass-dependent efficiency; 
is the cube of the radiative photon energy and reflects the expected energy dependence of the magnetic-dipole (Ml) 
matrix element; damping (E 7 ) describes a function to damp the diverging tail caused by the E^ dependence and is 
given in the form of E E()+ ^ E ,-e ) 2 as usec ^ ^y KEDR 22], where E is the peak energy of the transition photon. 

The experimental resolution function is determined from a signal MC sample with the width of the r\ c set to zero. A 
double Gaussian function is used for r\ c — > and a single Gaussian function for n c — > The mass-dependent 

efficiencies arc determined from phase-space MC samples. The background component in the channel n c — > is 
described by a third-order polynomial function. The background in the channel r/ c — > is composed of four parts: 

(1) contributions of J/ip — > S~S + , J/ip — > 7r°S - S + and J/ip — > £ A7r + 7r~ + c.c, with shapes and normalizations 
fixed in the fit; (2) a third-order Chebychev polynomial function representing the phase-space background contribution 
from J/ip — 7S~S + and other possible processes, with parameters set free in the fit. 

The signal detection efficiency is determined with MC simulated events by comparing the number of events after 
the event selection with the number of generated events. In the simulation, the decay J /ip — > 7?y c is generated using 



the helicity amplitude method 23j , and the radiative photon follows the angular distribution of 1 + cos 2 (8) , where 9 
is the polar angle of the radiative photon. The final state baryons' angular distributions are assumed to be uniformly 
distributed in the rest frame of the r\ c . 

The fitted curves are shown in Fig. [3] for r] c — > S + S _ and rj c — > S _ S + , where the mass and width of the r\ c are 
fixed to the newly measured results from BESIII [24 1 . A possible interference between the rj c resonance amplitude and 
the non- resonant background is neglected. The observed number of events, N b s , are listed in Table |U The statistical 
significances of the signals are calculated using the changes in the log-likelihood values and the number of degrees 
of freedom of the fits with and without the r\ c signal assumptions. For r/ c — s> the change in —ln(C) with 

A(d.o.f.) — 1 is 43.2, corresponding to a statistical significance of 9.3a. For r\ c — > the change in —ln(C) with 

A(d.o.f.) = 1 is 20.2, corresponding to a statistical significance of 6.4cr. The branching fraction of ?/ c — > is 
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calculated with: 



N , 



■bs 



peaking 



Nj^ x B(J/ip -> 777c) x /3 2 (S+ -> pvr°) x B 2 (ir° -> 77) x e' 
where N pea ki ng is the number of peaking background events determined from the background study, Nj/^ is the 



1J], B(J/^ -> 777c), -> p7T°) and 



total number of J ftp events, which is 2.25 x 10 8 with an uncertainty of 1.2% 
B(n° — > 77) are the branching fractions of J/ip — > jr) c , S + — > pn and 7r° — > 77, respectively [1|], and e is the total 
detection efficiency. The branching fraction of r\ c — > S~S + is calculated with: 



B( Vc -> S-S+) = 



AT. 



peaking 



N JN , x B(J/V -> 77? c ) x B 2 (E- -> Att") x /3 2 (A ^pvr-) x e' 
where /3(S~ — > Att - ) and £>(A p7r~) are the branching fractions of S~ — > and A — > pn~ , respectively 
results are summarized in Table HI 



1]. The 



TABLE I: Branching fractions of r\ c — » £ + E~ and ?7 C — >• H~H + obtained from this analysis and the predictions based on IML. 
For the measured branching fractions, the first uncertainty is statistical, the second experimental systematic, and the third is 
from input branching fractions taken from Ref. 



unc 

H 





7/c ->■ E+£" 


77c ->■ H"H + 


Statistical 


significance 


9.3(7 


6.4a 


N obs 




112 ± 15 


78 ± 14 


N peaking 




0.7 


2.0 


E 




5.3% 


5.5% 


Branching 


fraction (1CT 3 ) 


2. 11 ±0.28 ±0.18 ±0.50 


0.89 ±0.16 ±0.08 ±0.21 


Branching 


fraction based on IML [13] (1CT 3 ) 


0.51 - 1.00 


0.48 - 0.96 



VI. SYSTEMATIC UNCERTAINTIES 

The sources of systematic uncertainties for the two measurements are mainly from errors in the branching fractions 
of the known intermediate decay modes; the reconstruction and identification efficiencies of charged particles; the 
photon reconstruction; the ir°, E + , A and S - selection; vertex fits and kinematic fits; the fitting to the invariant-mass 
distributions; event generators and the total number of the J/tp events. The contributions are summarized in Table [TT1 

The tracking and identification efficiency of protons from the S + decay is determined using the J/ip — > Y, + Air~ data 
sample. The recoiling mass distribution of A7r~ pairs is fitted to obtain the E + signal yield, and the ratio between the 
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yields with and without the requirement of tracking and identifying the proton from the S + decay is determined. The 
tracking and PID efficiency for simulated MC events agrees within 2.0% with that obtained from the experimental 
data for each charged track. Hence, adding the uncertainties of the proton and anti-proton in quadrature, 2.8% is 
taken as the systematic uncertainty from reconstructing the final state charged tracks and their identification for 

The tracking and PID efficiencies of p, p, ir + and tt~ from S~ and S + decays are determined from analyzing 
J/iji — > — > AA7r + 7r~ — > pp7T + 7r + 7r~7r~ using a missing track method. Events are selected requiring all the 
tracks to be reconstructed except the one to be studied, and the invariant mass of the missing track predicted from 
the reconstructed tracks must be consistent with the invariant mass of the track to be studied. The tracking efficiency 
is then the fraction of the selected events with at least one additional track. The PID efficiency is obtained via the 
same missing track method. The tracking efficiency for MC simulated events is found to agree with that determined 
using data within 2.0% for each p, p track and 1.0% for each n + and tt~ track. Adding the uncertainties from p, p, 
7r + s and 7T~s in quadrature, 4.0% is taken as the systematic uncertainty for the six charged track final states. The PID 
efficiency for MC simulated events agrees with that determined using the data within 1.0% for each p, 2.0% for each 
p and 0.5% for each 7r + and 7r~, so 2.6% is taken as the systematic uncertainty for the pp7r + 7r + 7r~7r~ identification 
by adding the uncertainties in quadrature. 

The photon reconstruction efficiency is studied via three different methods: the missing photon method, the missing 
7T° method and the 7r° decay angle method with ip' — > ir + Tr~J/ijj — > ir + ir~ p°7r° , ip' — > tt tt q J/i(j — s> tt q tt°1 + 1" and 
J/ip — >■ p°tt° events, respectively. The efficiency difference between data and MC simulated events is within 1.0% 



for each photon 25j. Thus, 5.0% and 1.0% are taken as the systematic uncertainty due to photon reconstruction for 
i] c — > and i] c — > whose final states contain five photons and one photon, respectively. 

The uncertainty of the ir° selection is determined with the data sample J ftp — > E~A7r + — > Tr a ppir + ir~ . The w° 
selection efficiency is determined from the change in the S~ signal yield from fitting the An + recoiling mass distribution 
with and without the ir° selection requirement. The difference between beam data and MC simulated events on the 
7r°-selection efficiency is within 0.5% per 7r°; hence 1.0% is taken as the systematic uncertainty from w° selection for 
7? c ->£+£-. 

Samples of J/ij) -} jK* + K*~ — > ~fK + K~ir a ir°, J/ip — > ppr/ — > pp7r 7r°7r° and J/ip — > jr) c — > jK + K~7r + ir + ir~ir~ 
arc selected to study the efficiency difference between beam data and simulated MC events in the kinematic fitting 
analysis for r\ c -> £+£~ and rj c ->• S~S+. In r) c -> £+£~, the sample of J/ip — > r )K* + K*~ — > ~fK + K~Tr°Tr is selected 
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to estimate the efficiency of the first two xlc requirements: X4c(PP 7r ° 7r °) > 200 and xlci^PP^ ' 7r °) < X4c(^PP 7r ° 7r °)' 
and the efficiency of the xXc^CIPP 1 ^ ' 7r °) < 30 requirement is estimated by the change in the r\ signal yield from fitting 
the pp recoiling mass distribution from J/tp — > ppr\ — > pp7r 7r°7r° when the \\c °^ the •W' ~~ * PP"" "" "" hypothesis 
is less than 30. In 77 c — > we select a clean J/ip — > -fr] c — > jK + K~ir + ir + ir~7r~ sample, plot the 4C kinematic 

fitting efficiency at different \\ c requirements and obtain the efficiency for the requirements as described in the event 
selection section. The estimated systematic uncertainties are 4.3% and 3.8% from kinematic fitting for rj c — > E + E~ 
and rj c — > S~S + , respectively. 

The uncertainty from the E + -mass window requirement is estimated by selecting a sample of J/ip — > E + S _ events 
and by studying the efficiency difference between beam data and simulated MC events. An uncertainty of 0.6% is 
found. 

The uncertainties from the vertex fits and from the S~ , A-mass window requirements are estimated from a sample 
of J/ip — > — > AA7r + 7r~ — > ppir + TT + Tr~Tr~ events. The efficiency difference between beam data and simulated 
MC events is within 0.6%, 0.3% and 0.3% for the vertex fits, S~ and A-mass window requirements, respectively. 

Uncertainties from event generators are studied by comparing results with different models that were used for the 
generation of the signal events. The decays r\ c — > £ + E~ and r\ c — > are generated with another model using 

the helicity amplitude, and assuming that the baryons are uniformly distributed in the rest frame of ?y c ; the decays 
E + — > p7r°, S~ — > Att~ and A — > pir~ are generated with another model, which takes parity violation effects into 
consideration. The efficiency differences are 0.4% and 2.8% for r\ c — > and r\ c — > respectively. 

Uncertainties from fitting the invariant-mass distributions of E + S _ and S _ S + pairs are estimated by varying signal 
md ,» — „« pes a»d 1. cor— fitting mg e. Tll e ™ 8 and ^ of th e , k „ rarM by u a«o rfills 

to the new measurements from BESIII [24|; the damping function is changed from the form used by KEDR 

- ~" I 
e with (3 fixed at 65 MeV, which was used by CLEO [26J; the MC signal shape is convoluted with a Gaussian 

with the width as a free parameter in the fit to study a possible uncertainty from the mass resolution determined 

from simulated MC events; the background shapes are varied either through the order of the polynomial or the 

normalization of fixed parts; the fitting range is varied to either a narrower or a wider one. Taking all the factors 

described above into account and by adding the uncertainties from each factor in quadrature, the uncertainties due 

to the fitting procedures are estimated to be 4.7% and 6.4% for r\ c — > E + E _ and j] c — > S~S + , respectively. 

The measured branching fractions of the peaking background channels have uncertainties around ~ 20 — 30%. The 



m to 
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uncertainties from the number of peaking background events are estimated by assigning conservative estimates of 50% 
to the uncertainties of the measured branching fractions of r\ c — > ppTr°ir°, r\ c — > pp7r + 7r + 7r~7r~ and r\ c — > AA7r + 7r~. 

The total number of J/tjj events is determined from analyzing J ftp inclusive hadronic decays, and the uncertainty 
is 1.2% [l4|. 



Limited knowledge of the branching fractions, B(J/ip — > JTjc), £>(£+ — > P^°), and B(A — > pir~) contribute 23.5%, 
0.6%, and 0.8% uncertainty to the measured branching fractions, respectively 1|. The first of these is the dominant 
source of systematic uncertainty, as indicated in Table [XXJ 

All the systematic uncertainties and their sources for the channels rj c — > and rj c — > are summarized 

in Tabic HU The quadratic sum of all the systematic uncertainties that solely stem from our experiment are 8.7% and 
9.5% in the branching fraction measurements of r\ c — > and r\ c — > S~S + , respectively. The total systematic 

uncertainty is about 25% for both measurements. 



TABLE II: Systematic uncertainties (%) in the branching fraction measurements of rj c — >• E + E and rj, 



c —r — — 



Source 


r? c -> E+E^ 




Tracking and PID 


2.8 


4.8 


Photon reconstruction 


5.0 


1.0 


7T° selection 


1.0 




E + mass window 


0.6 




A mass window 




0.3 


H - mass window 




0.3 


Vertex fits 




0.5 


Kinematic fits 


4.3 


3.8 


Signal fitting 


4.7 


6.4 


Event generators 


0.4 


2.8 


Peaking background 


0.3 


1.3 




1.2 


1.2 


Intermediate states 


23.5 


23.6 


Total (BESIII) 


8.7 


9.5 


Total 


25.1 


25.5 
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VII. SUMMARY 



Using 2.25 x 10 8 J ftp events collected with the BESIII detector, the decays J/i/j — > ^r\ c — > 7£ + £ and J ftp —> 
jTjc — > 7S~S + are observed for the first time, and their branching fractions are measured to be: 

B(J/ip -> 7?7 C -> 7S+E-) = (3.60 ± 0.48 ± 0.31) x 10~ 5 , 
B(J/ip -> 7?/c -> 7S-S+) = (1.51 ± 0.27 ± 0.14) x 10- 5 . 



Using the known value of B(J/ip — > jrjc) = (1-7 ± 0.4)% 
are obtained: 



1], the branching fractions of r\ c —> E + S and r\ c — ► S s + 



S(?y c -> £+£") = (2.11 ± 0.28 ± 0.18 ± 0.50) x 10" 3 , 
B(r] c -> S"S+) = (0.89 ± 0.16 ± 0.08 ± 0.21) x 10" 3 , 

where the first uncertainties are statistical, the second systematic, and the third uncertainties are from the precision 
of the intermediate branching fractions. 



Table U compares the results of our measurements with the predictions from charmed- meson loop calculations 13 1 . 
The measured branching fraction of J] c — > is larger than the prediction, while the measured branching fraction 

of rj c — > 5 _ S + agrees with the prediction. Among the four r\ c baryonic decays (rj c — > pp, A A, £ + S~, and S - EE + ), 
only rj c — > disagrees with the prediction, which may indicate the violation of SU(3) symmetry. 

The precision of the branching fraction measurements of rj c — > and rj c — > are limited by statistics, and 

the dominating systematic error stems from the uncertainty in the branching fraction of J/ip — > "yrjc, which cannot 
be reduced without a thorough theoretical understanding of the r\ c line shape in Ml transitions in the charmonium 
system. 
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